Measurements of single top quark production cross sections and | Vtt, | in pp collisions 

at y/s = 1.96 TeV 



V.M. Abazov, 34 B. Abbott, 72 B.S. Acharya, 28 M. Adams, 48 T. Adams, 46 G.D. Alexeev, 34 G. Alkhazov, 38 
A. Alton , 60 G. Alverson, 59 GA. Alves, 2 M. Aoki, 47 M. Arov, 57 A. Askew, 46 B. Asman, 40 S. Atkins, 57 

0. Atramentov, 64 K. Augsten, 9 C. Avila, 7 J. BackusMayes, 79 F. Badaud, 12 L. Bagby 47 B. Baldin, 47 
D.V. Bandurin, 46 S. Banerjee, 28 E. Barberis, 59 P. Baringer, 55 J. Barreto, 3 J.F. Bartlett 47 U. Bassler, 17 
V. Baztcrra 48 A. Bean, 55 M. Begalli, 3 M. Begel, 70 C. Belanger-Champagnc 40 L. Bcllantoni, 47 S.B. Beri, 26 

G. Bcrnardi, 16 R. Bernhard, 21 I. Bertram, 41 M. Besangon, 17 R. Beuselinck 42 V.A. Bezzubov, 37 RC. Bhat, 47 
V. Bhatnagar, 26 G. Blazey 49 S. Blessing 46 K. Bloom, 63 A. Boehnlein, 47 D. Boline, 69 E.E. Boos, 36 G. Borissov, 41 

T. Bose, 58 A. Brandt, 75 O. Brandt, 22 R. Brock, 61 G. Brooijmans, 67 A. Bross 47 D. Brown, 16 J. Brown, 16 X.B. Bu, 47 
M. Buehler 47 V. Buescher, 23 V. Bunichev, 36 S. Burdin 6 , 41 T.H. Burnett, 79 CP. Buszello 40 B. Calpas, 14 
E. Camacho-Perez, 31 M.A. Carrasco-Lizarraga, 55 B.CK. Casey 47 H. Castilla-Valdez, 31 S. Chakrabarti, 69 
D. Chakraborty 49 K.M. Chan, 53 A. Chandra, 77 E. Chapon, 17 G. Chen, 55 S. Chevalier-Thery, 17 D.K. Cho, 74 
S.W. Cho, 30 S. Choi, 30 B. Choudhary, 27 S. Cihangir, 47 D. Claes, 63 J. Clutter, 55 M. Cooke 47 W.E. Cooper 47 
M. Corcoran, 77 F. Coudcrc, 17 M.-C. Cousinou, 14 A. Croc, 17 D. Cutts, 74 A. Das, 44 G. Davics 42 K. De, 75 
S.J. dc Jong, 33 E. De La Cruz-Burclo, 31 F. Deliot, 17 M. Demarteau, 47 R. Demina, 68 D. Denisov 47 S.P. Dcnisov, 37 
S. Dcsai, 47 C. Deterre, 17 K. DeVaughan, 63 H.T. Diehl 47 M. Diesburg 47 P.F. Ding 43 A. Dominguez, 63 T. Dorland, 79 
A. Dubey, 27 L.V. Dudko, 36 D. Duggan, 64 A. Duperrin, 14 S. Dutt, 26 A. Dyshkant, 49 M. Eads, 63 D. Edmunds, 61 
P. Ellcr' 1 48 J. Ellison, 45 V.D. Elvira 47 Y. Enari, 16 H. Evans, 51 A. Evdokimov, 70 V.N. Evdokimov, 37 G. Facini, 59 
T. Fcrbel, 68 F. Fiedler, 23 F. Filthaut, 33 W. Fisher, 61 H.E. Fisk 47 C. Fockc' 1 48 M. Fortncr, 49 H. Fox 41 S. Fuess 47 
A. Garcia-Bellido, 68 G.A Garci'a-Guerra c , 31 V. Gavrilov, 35 P. Gay, 12 W. Ceng, 14 - 61 D. Gerbaudo, 65 C.E. Gerber, 48 
Y. Gcrshtcin, 64 G. Ginther , 47 - 68 G. Golovanov, 34 A. Goussiou, 79 P.D. Grannis, 69 S. Grcdcr, 18 H. Greenlee, 47 
Z.D. Greenwood, 57 E.M. Gregores, 4 G. Grenier, 19 Ph. Gris, 12 J.-F. Grivaz, 15 A. Grohsjcan, 17 S. Griinendahl, 47 
M.W. Griinewald, 29 T. Guillemin, 15 G. Gutierrez 47 P. Gutierrez, 72 A. Haas d , 67 S. Hagopian, 46 J. Haley, 59 
L. Han, 6 K. Harder, 43 A. Harcl, 68 J.M. Hauptman, 54 J. Hays, 42 T. Head, 43 T. Hebbekcr, 20 D. Hedin, 49 

H. Hegab, 73 A.P. Heinson 45 U. Heintz, 74 C. Hensel, 22 I. Heredia-De La Cruz, 31 K. Herner, 60 G. Hesketh 6 43 
M.D. Hildreth, 53 R. Hirosky, 78 T. Hoang, 46 J.D. Hobbs, 69 B. Hoeneisen, 11 M. Hohlfeld, 23 Z. Hubacek, 9 - 17 
N. Huske, 16 V. Hynek, 9 I. Iashvili, 66 Y. Ilchenko, 76 R. Illingworth, 47 A.S. Ito 47 S. Jabeen, 74 M. Jaffre, 15 

D. Jamin, 14 A. Jayasinghe, 72 R. Jesik, 42 K. Johns, 44 M. Johnson, 47 A. Jonckheere, 47 P. Jonsson, 42 J. Joshi, 26 
A.W. Jung 47 A. Juste, 39 K. Kaadze, 56 E. Kajfasz, 14 D. Karmanov, 36 P.A. Kasper 47 I. Katsanos, 63 R. Kehoe, 76 
S. Kermiche, 14 N. Khalatyan, 47 A. Khanov, 73 A. Kharchilava, 66 Y.N. Kharzheev, 34 J.M. Kohli, 26 A.V. Kozelov, 37 

J. Kraus, 61 S. Kulikov, 37 A. Kumar, 66 A. Kupco, 10 T. Kurca, 19 V.A. Kuzmin, 36 J. Kvita, 8 S. Lammers, 51 
G. Landsberg, 74 P. Lebrun, 19 H.S. Lee, 30 S.W. Lee, 54 W.M. Lee, 47 J. Lellouch, 16 L. Li 45 Q.Z. Li 47 S.M. Lietti, 5 
J.K. Lim, 30 D. Lincoln 47 J. Linnemann, 61 V.V. Lipaev, 37 R. Lipton, 47 Y. Liu, 6 A. Lobodenko, 38 M. Lokajicek, 10 
R. Lopes de Sa, 69 H.J. Lubatti, 79 R. Luna-Garcia^, 31 A.L. Lyon, 47 A.K.A. Maciel, 2 D. Mackin, 77 R. Madar, 17 

R. Magaha-Villalba, 31 S. Malik, 63 V.L. Malyshev, 34 Y. Maravin, 56 J. Martinez-Ortega, 31 R. McCarthy, 69 
C.L. McGivern, 55 M.M. Meijer, 33 A. Melnitchouk, 62 D. Menezes 49 P.G. Mercadante, 4 M. Merkin, 36 A. Meyer, 20 
J. Meyer, 22 F. Miconi, 18 N.K. Mondal, 28 G.S. Muanza, 14 M. Mulhearn, 78 E. Nagy, 14 M. Naimuddin, 27 M. Narain, 74 
R. Nayyar, 27 H.A. Ncal, 60 J. P. Negret, 7 P. Neustroev, 38 S.F. Novacs, 5 T. Nunnemann, 24 G. Obrant*, 38 J. Orduna, 77 
N. Osman, 14 J. Osta, 53 G.J. Otero y Garzon, 1 M. Padilla, 45 A. Pal, 75 N. Parashar, 52 V. Parihar, 74 S.K. Park, 30 
J. Parsons, 67 R. Partridge 6 *, 74 N. Parua, 51 A. Patwa, 70 B. Penning, 47 M. Perfilov, 36 K. Peters, 43 Y. Peters, 43 
K. Pctridis 43 G. Petrillo, 68 P. Pctroff, 15 R. Piegaia, 1 M.-A. Plcicr, 70 P.L.M. Podcsta-Lcrmaf, 31 V.M. Podstavkov, 47 
P. Polozov, 35 A.V. Popov, 37 M. Prcwitt, 77 D. Price, 51 N. Prokopcnko, 37 H.B. Prosper, 46 S. Protopopescu, 70 
J. Qian, 60 A. Quadt, 22 B. Quinn, 62 M.S. Rangel, 2 K. Ranjan, 27 P.N. Ratoff, 41 I. Razumov, 37 P. Renkel, 76 
M. Rijssenbeek, 69 I. Ripp-Baudot, 18 F. Rizatdinova, 73 M. Rominsky 47 A. Ross, 41 C. Royon, 17 P. Rubinov 47 
R. Ruchti, 53 G. Safronov, 35 G. Sajot, 13 P. Salcido 49 A. Sanchez-Hernandez, 31 M.P. Sanders, 24 B. Sanghi 47 
A.S. Santos, 5 G. Savage, 47 L. Sawyer, 57 T. Scanlon 42 R.D. Schamberger, 69 Y. Scheglov, 38 H. Schellman, 50 
T. Schliephakc, 25 S. Schlobohm, 79 C. Schwanenberger, 43 R. Schwienhorst, 61 J. Sekaric, 55 H. Severini, 72 
E. Shabalina, 22 V. Shary, 17 A. A. Shchukin, 37 R.K. Shivpuri, 27 V. Simak, 9 V. Sirotenko, 47 P. Skubic, 72 P. Slattery, 68 
D. Smirnov, 53 K.J. Smith, 66 G.R. Snow, 63 J. Snow, 71 S. Snyder, 70 S. Soldner-Rembold, 43 L. Sonnenschein, 20 
K. Soustruznik, 8 J. Stark, 13 V. Stolin, 35 D.A. Stoyanova, 37 M. Strauss, 72 D. Strom, 48 L. Stuttc 47 L. Suter, 43 
P. Svoisky, 72 M. Takahashi, 43 A. Tanasijczuk, 1 M. Titov, 17 V.V. Tokmenin, 34 N. Triplett, 54 Y.-T. Tsai, 68 



2 



K. Tschann-Grimm, 69 D. Tsybychev, 69 B. Tuchming, 17 C. Tully, 65 L. Uvarov, 38 S. Uvarov, 38 S. Uzunyan, 49 
R. Van Kooten, 51 W.M. van Lccuwen, 32 N. Varelas, 48 E.W. Varncs, 44 LA. Vasilyev, 37 P. Verdier, 19 
L.S. Vcrtogradov, 34 M. Vcrzocchi, 47 M. Vesterinen, 43 D. Vilanova, 17 P. Vokac, 9 H.D. Wahl, 46 M.H.L.S. Wang, 47 

J. Warchol, 53 G. Watts, 79 M. Wayne, 53 M. Weber' 1 , 47 L. Wclty-Rieger, 50 A. White, 75 D. Wicke, 25 
M.R.J. Williams, 41 G.W. Wilson, 55 M. Wobisch, 57 D.R. Wood, 59 T.R. Wyatt, 43 Y. Xie, 47 C. Xu, 60 S. Yacoob, 50 
R. Yamada, 47 W.-C. Yang, 43 T. Yasuda, 47 Y.A. Yatsunenko, 34 Z. Ye, 47 H. Yin, 47 K. Yip, 70 S.W. Youn 47 
J. Yu, 75 S. Zclitch, 78 T. Zhao, 79 B. Zhou, 60 J. Zhu, 60 M. Zielinski, 68 D. Zieminska, 51 and L. Zivkovic 74 

(The DO Collaboration*) 

I Umversidad de Buenos Aires, Buenos Aires, Argentina 
2 LAFEX, Centra Brasileiro de Pesquisas Fisicas, Rio de Janeiro, Brazil 

3 Universidade do Estado do Rio de Janeiro, Rio de Janeiro, Brazil 
4 Universidade Federal do ABC, Santo Andre, Brazil 
5 Instituto de Fisica Teorica, Universidade Estadual Paulista, Sao Paulo, Brazil 
6 University of Science and Technology of China, Hefei, People's Republic of China 
Universidad de los Andes, Bogota, Colombia 
8 Charles University, Faculty of Mathematics and Physics, 
Center for Particle Physics, Prague, Czech Republic 
9 Czech Technical University in Prague, Prague, Czech Republic 
10 Center for Particle Physics, Institute of Physics, 
Academy of Sciences of the Czech Republic, Prague, Czech Republic 

II Universidad San Francisco de Quito, Quito, Ecuador 
12 LPC, Universite Blaise Pascal, CNRS/IN2P3, Clermont, France 

1S LPSC, Universite Joseph Fourier Grenoble 1, CNRS/IN2P3, 
Institut National Polytechnique de Grenoble, Grenoble, France 
U CPPM, Aix-Marseille Universite, CNRS/IN2P3, Marseille, France 
15 LAL, Universite Paris-Sud, CNRS/IN2P3, Orsay, France 
ie LPNHE, Unwersites Paris VI and VII, CNRS/IN2P3, Paris, France 
"CEA, Irfu, SPP, Saclay, France 
18 IPHC, Universite de Strasbourg, CNRS/IN2P3, Strasbourg, France 
19 IPNL, Universite Lyon 1, CNRS/IN2P3, Villeurbanne, France and Universite de Lyon, Lyon, France 
20 III. Physikalisches Institut A, RWTH Aachen University, Aachen, Germany 
21 Physikalisches Institut, Universitdt Freiburg, Freiburg, Germany 
22 II. Physikalisches Institut, Georg- August- Universitdt Gdttmgen, Gdttmgen, Germany 
2S Institut fur Physik, Universitdt Mainz, Mainz, Germany 
2Jl Ludwig-Maximilians-Universitat Munchen, Miinchen, Germany 
25 Fachbereich Physik, Bergische Universitdt Wuppertal, Wuppertal, Germany 
26 Panjab University, Chandigarh, India 
21 Delhi University, Delhi, India 
28 Tata Institute of Fundamental Research, Mumbai, India 

29 University College Dublin, Dublin, Ireland 
Korea Detector Laboratory, Korea University, Seoul, Korea 
31 CINVESTAV, Mexico City, Mexico 
J2 Nikhef, Science Park, Amsterdam, the Netherlands 
33 Radboud University Nijmegen, Nijmegen, the Netherlands and Nikhef, Science Park, Amsterdam, the Netherlands 

3Jt Joint Institute for Nuclear Research, Dubna, Russia 
35 Institute for Theoretical and Experimental Physics, Moscow, Russia 
36 Moscow State University, Moscow, Russia 
37 Institute for High Energy Physics, Protvino, Russia 
38 Petersburg Nuclear Physics Institute, St. Petersburg, Russia 
39 Institucio Catalana de Recerca i Estudis Avangats (ICREA) and Institut de Fisica d'Altes Energies (IFAE), Barcelona, Spam 
40 Stockholm University, Stockholm and Uppsala University, Uppsala, Sweden 
J>1 Lancaster University, Lancaster LAI 4YB, United Kingdom 
^Imperial College London, London SW7 2AZ, United Kingdom 
43 The University of Manchester, Manchester MIS 9PL, United Kingdom 

JfJt University of Arizona, Tucson, Arizona 85721, USA 
45 University of California Riverside, Riverside, California 92521, USA 

40 Florida State University, Tallahassee, Florida 32306, USA 
47 Fermi National Accelerator Laboratory, Batavia, Illinois 60510, USA 
48 University of Illinois at Chicago, Chicago, Illinois 60607, USA 
49 Northern Illinois University, DeKalb, Illinois 60115, USA 
50 Northwestern University, Evanston, Illinois 60208, USA 
51 Indiana University, Bloommgton, Indiana 47405, USA 



3 



52 Purdue University Calumet, Hammond, Indiana 46323, USA 
53 University of Notre Dame, Notre Dame, Indiana 46556, USA 
5J, Iowa State University, Ames, Iowa 50011, USA 
55 University of Kansas, Lawrence, Kansas 66045, USA 
56 Kansas State University, Manhattan, Kansas 66506, USA 
5 Louisiana Tech University, Ruston, Louisiana 71272, USA 

58 Boston University, Boston, Massachusetts 02215, USA 
59 Northeastern University, Boston, Massachusetts 02115, USA 
60 University of Michigan, Ann Arbor, Michigan 48109, USA 
61 Michigan State University, East Lansing, Michigan 48824, USA 
62 University of Mississippi, University, Mississippi 38677, USA 
63 University of Nebraska, Lincoln, Nebraska 68588, USA 
b ' ''Rutgers University, Piscataway, New Jersey 08855, USA 
65 Princeton University, Princeton, New Jersey 08544, USA 
66 State University of New York, Buffalo, New York 14260, USA 
61 Columbia University, New York, New York 10027, USA 
68 University of Rochester, Rochester, New York 14627, USA 
69 State University of New York, Stony Brook, New York 11794, USA 
70 Brookhaven National Laboratory, Upton, New York 11973, USA 
1 Langston University, Langston, Oklahoma 73050, USA 
72 University of Oklahoma, Norman, Oklahoma 73019, USA 
73 Oklahoma State University, Stillwater, Oklahoma 74078, USA 
7J < Brown University, Providence, Rhode Island 02912, USA 

75 University of Texas, Arlington, Texas 76019, USA 
7S Southern Methodist University, Dallas, Texas 75275, USA 
Rice University, Houston, Texas 77005, USA 

78 University of Virginia, Charlottesville, Virginia 22901, USA 

79 University of Washington, Seattle, Washington 98195, USA 

(Dated: August 15, 2011) 

We present measurements of production cross sections of single top quarks in pp collisions at y/s = 
1.96 TeV in a data sample corresponding to an integrated luminosity of 5.4 fb _1 collected by the 
DO detector at the Fermilab Tevatron Collider. We select events with an isolated electron or muon, 
an imbalance in transverse energy, and two, three, or four jets, with one or two of them containing 
a bottom hadron. We obtain an inclusive cross section of o(pp — > tb + X, tqb + X) = 3.43 ±0.74 pb 
and use it to extract the CKM matrix element 0.79 < |Vti,| < 1 at the 95% C.L. We also measure 
a(pp — ¥ tb + X) — 0.68 ±o;l| pb and o(pp — ¥ tqb + X) = 2.86±o;e3 pb when assuming, respectively, 
tqb and tb production rates as predicted by the standard model. 

PACS numbers: 14.65.Ha; 12.15.Ji; 13.85.Qk; 12.15.Hh 



I. INTRODUCTION 



Top quarks are produced at hadron colliders as tt pairs 
via the strong interaction or singly via the electroweak 
interaction [J 0. Because of the larger production rate 
and higher signal-to-background ratio, the production 
of tt pairs is better studied and indeed it was through 
the tt production process that the existence of the top 
quark was established in 1995 at the Fermilab Tevatron 
Collider 0, 0] . The observation of the single top quark 
production, however, was possible after CDF and DO 
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collaborations accumulated « 50 times more integrated 
luminosity than what was needed for observation of top 
quarks in tt production 0, @] ■ Single top quark events are 
produced at about half of the rate of top quark pairs and 
with lower jet multiplicities, and therefore their study 
is more susceptible to contamination from background 
processes. 

Electroweak production of top quarks at the Tevatron 
proceeds mainly via the decay of a time-like virtual W 
boson accompanied by a bottom quark in the s channel 
(tb = tb+tb) Q or via the exchange of a space- like virtual 
W boson between a light quark and a bottom quark in 
the t channel (tqb = tqb + tqb, where q refers to the light 
quark or antiquark) Figure[T]shows the lowest level 

Feynman diagrams for s- and i-channel production [Ioj |. 
A third process tW, in which the top quark is produced 
together with a W boson, has a small cross section at 
the Tevatron Q and is therefore not considered in this 
analysis. 

Single top quark events can be used to probe the 
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FIG. 1: [color online] Lowest level Feynman diagrams for (a) 
tb and (b) tqb single top quark production. 



Wtb vertex and to directly measure the mag nitude 
of the Cabibbo-Kobayashi-Maskawa (CKM) [ljj quark 
mixing matrix element \Vtb\- Under the assumptions 
that there are only three quark generations and that 
the CKM matrix is unitary, the matrix elements are 
severely constrained QJ]: \V td \ = (8.62±g;|g) x lfT 3 , 



(4.0318$) x 10- 



and \V t 



lb 



0.999152 



+0.000030 
-0.000045- 



\Vts\ 

However, in several extensions of the standard model 
(SM) involving, for instance, a fourth generation of 
quarks or an additional heavy quark singlet that mixes 
with the top quark, \Vtb\ can be significantly smaller 
than unity |l3[ . A direct determination of | Vtb | , without 
assuming unitarity or three generations, is possible 
through the measurement of the total single top quark 
production cross section [l4j . The current measured 
value for the total single top quark cross section is 



+0.58 
-0.47 



pb at yfs = 1.96 TeV, resulting in |Vtj| 



2.76 

0.88 ± 6.07 with a limit of \V tb \ > 0.77 at the 95% C.L. 
assuming a top quark mass m t = 170 GeV (l5j . 

Previous measurements of single top quark production 
cross sections included events from both the 

tb and tqb processes, assuming a ratio of cross sections Q 
for the two processes based on the SM. However, several 
beyond-the-SM theories predict individual tb and tqb 
cross sections that deviate from the SM. Examples 
include models with additional quark generations fl3| . 
new heavy bosons [20| . flavor-changing neutral currents 
(FCNC) [2lj ]. or anomalous top quark couplings [22T - |24j ] . 
It is therefore important to also measure the individual 
tb and tqb production rates. 

Using data corresponding to 5.4 fb _1 of integrated 
luminosity recorded with the DO detector [25| , we present 
an improved measurement of the production rate of 
tb+tqb. We also present measurements of the production 
rates of the individual tb and tqb processes performed 
assuming, respectively, tqb and tb production rates as 
predicted by the SM. Finally, we present a new direct 
measurement of | Vtb I extracted from the measured tb+tqb 
cross section. 



II. EVENT SELECTION 

This analysis extends previous work by the DO Collab- 
oration 0, O, [It], HH and uses the same data, event 
selection, and modeling of signal and background as in 
Ref. [I?]] ■ It differs however, in the assumptions used 
to extract the cross sections of the individual tb and tqb 
production modes. 

The data were collected with a logical OR of many 
trigger conditions, which together are fully efficient for 
the single top quark signal. We select events containing 
only one isolated electron or muon with high trans- 
verse momentum (pr) and having a large imbalance in 
the transverse energy ($t) indicative of the presence of 
a neutrino. Events originating from single top quark 
production are expected to contain at least one b quark 
jet from the decay of the top quark and a second b quark 
jet in the s channel, or a light quark jet and a spectator b 
quark jet for the t channel. In both cases, gluon radiation 
can give rise to additional jets. Events are selected 
with two, three, or four jets reconstructed using a cone 
algorithm [28[ in (y, (j>) space, where y is the rapidity 
and 4> is the azimuthal angle, and the cone radius 0.5. 
The jets must satisfy the following conditions: leading 
jet px > 25 GeV, other jets with px > 15 GeV, and with 
pseudorapidities of all jets \r]\ < 3.4. Requirements are 
also placed on Jfi T ■ 20 < $t < 200 GeV for events with 
two jets, and 25 < $t < 200 GeV for events with three or 
four jets. The maximum $t requirement removes events 
that suffer from poor modeling of the high energy tail of 
the muon momentum resolution. We require one isolated 
electron with |?y| < 1.1 and px > 15 (20) GeV for events 
with two (three or four) jets, or one isolated muon with 
\r]\ < 2.0 and p T > 15 GeV. 

The sample resulting from this selection is dominated 
by W bosons produced in association with jets (VF+jets), 
with smaller contributions from ti pairs decaying into 
the single lepton plus jets final state or the dilepton 
final state when one lepton or some jets are not 
reconstructed. Multijet events also contribute to the 
background when a jet is misidentified as an isolated 
electron or a heavy-flavor quark decays to a muon that 
satisfies isolation criteria, in combination with misre- 
construction of # T . Diboson (WW, WZ, ZZ) and 
2+jets processes contribute only marginally to the total 
background but are taken into account. The background 
from multijets is kept small (< 6%) by requiring that the 
total scalar sum (i?r) of the transverse momenta of the 
final-state objects (lepton, and jets) be larger than 
120/140/160 GeV for events with 2 /3 /4 jets, and that the 
JpT does not point along the transverse direction of the 
lepton or the leading jet. Soft-scattering processes are 
suppressed by requiring a minimum value for the total 
scalar sum of the px of the jets [Ht (alljets)] ranging from 
50 to 100 GeV, depending on the number of jets in an 
event and the data collection period. 

To enhance the signal fraction, one or two of the jets 
are required to contain long lived bottom hadrons (b 
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jets), as determined through a multivariate 6-tagging 
algorithm [30] . This algorithm uses several variables to 
discriminate 6 jets from other jets such as: (i) decay 
length significance of the secondary vertex, (ii) the x 2 
per degree of freedom of the secondary vertex fit, (iii) 
weighted combination of the tracks' impact parameter 
significances, (iv) probability that the jet originates from 
the primary pp interaction vertex, (v) number of tracks 
used to reconstruct the secondary vertex, (vi) mass of 
the secondary vertex, and (vii) number of secondary 
vertices found inside the jet. To improve sensitivity 
to signal, the samples are divided into six independent 
analysis channels, depending on the jet multiplicity (two, 
three, or four jets), and the number of 6-tagged jets 
(one or two). The efficiency of the event selection, 
including branching fraction and the 6-tagging require- 
ments, is (2.9 ±0.4)% for tb and (2.0 ±0.3)% for tqb. The 
tqb process has a lower acceptance than the tb channel 
because the second 6-jet has low transverse momentum 
and is difficult to identify. We apply additional require- 
ments to select two control samples used to test whether 
the background model reproduces the data in regions 
dominated by one specific type of background. The 
control sample dominated by W±jets is required to have 
exactly two jets, Ht < 175 GeV, and only one 6-tagged 
jet where VF±jets events constitute 82% of this sample, 
and the ti component is less that 2%. The control sample 
dominated by ti is required to have exactly four jets, 
Ht > 300 GeV, and one or two 6-tagged jets where ti 
events constitute 84% of the sample, and the Ty±jets 
component is 12%. 



III. MODELS FOR SIGNAL AND 
BACKGROUND 

Single top quark events are modeled for a top quark 
mass nit = 172.5 GeV using the COMPHEP-based effective 
next-to-leading order (NLO) Monte Carlo (MC) event 
generator SINGLETOP [31j . which preserves spin infor- 
mation in the decays of the top quark and the W boson 
and provides event kinematics that reproduce distri- 
butions predicted by NLO calculations [H, [H[. The 
ti, W+jets, and Z+jets events are simulated with the 
ALPGEN leading orden MC generator [34j|. Diboson 
processes are modeled using pythia (35|. For all these 
MC samples, pythia is also used to evolve parton 
showers and to model proton remnants and hadronization 
of all generated partons. The presence of additional pp 
interactions is modeled by events selected from random 
beam crossings matching the instantaneous luminosity 
profile in the data. All MC events are passed through a 
GEANT-based simulation [36| of the DO detector. 

Differences between simulation and data in lepton 
and jet reconstruction efficiencies and resolutions, jet 
energy scale, and 6-tagging efficiencies are corrected in 
the simulation by applying correction functions measured 
from separate data samples. Comparisons of ALPGEN 



with data and with other generators show small discrep- 
ancies in distributions of jet pseudorapidity and angular 
separations between jets [331 . We therefore correct 
the ALPGEN PF±jets and Z+jets samples by sequen- 
tially applying polynomial reweighting functions param- 
eterized by the leading and second-leading jet r\, AR = 
y/ (A</>) 2 + (A?/) 2 between the two leading jets, and third- 
and fourth-leading jet w, if applicable. These functions 
are derived from the ratio between the number of W+jets 
and Z+jets events observed in data and the event 
yields predicted by MC. After these corrections, the MC 
description is in good agreement with our high statistics 
sample of events prior to the application of 6-tagging. 
The multijet background is modeled using the selection 
discussed in SecHIl but choosing events that fail isolation 
criteria for leptons. 

MC samples are scaled to the theoretical cross section 
at approximately NNLO [l| for ti, and NLO [HI for 
Z+jets and diboson cases. The contributions from 
tV~+jets and multijet are normalized by comparing the 
prediction for background to data before 6-tagging. We 
use a procedure that relies on three distributions [lepton 
Pt, $t and W reconstructed mass in the transverse plane 
Mt{W)] that have distinctive shapes for VF+jets and 
multijets events and are thus sensitive to their relative 
contributions in the selected sample. The normalization 
scale factors for W+jets (Awjets) and multijet (A mu itij e ts) 
are constrained by the following equation: 

N = AyVjets^Wjets + Amultijets-^multijets, (1) 
Where N = A da ta - Nnon-Wjcts and TVdata, Nnon-Wjcts, 

■^Wjets, and iV mu itijets are the event yields in data, non- 
FF+jets MC, VF±jets, and multijet samples, respectively. 
The FF+jets sample contains events with light flavor 
(Wjj, j — u, d, s) and heavy flavor (Wjc, Wcc, and Wbb) 
quarks. The non-14 7 ±jets MC samples include single 
top quark, ti, Z+jets, and diboson production. The 
values of Awjets and A mu itij e ts are varied to maximize the 
product of the Kolmogorov-Smirnov test values [39[ for 
the three kinematic distributions. This procedure is done 
separately for events with two, three, and four jets and 
for each lepton flavor. After the normalization, the total 
sum of the TF±jets and multijets yields plus the small 
contributions from ti, single top, Z+jets, and diboson 
production equals the total data yield for each of the six 
analysis channels. 

Without modifying the overall normalization of the 
W+jets MC sample, we apply an additional scale factor 
to W and Z boson events produced in conjunction with 
heavy-flavor jets (6 or c) to match NLO calculations [38| : 
Wbb and Wcc by 1.47, Zbb by 1.52, Zee by 1.67, and 
Wcj by 1.32. We evaluate whether an additional normal- 
ization factor Xhf is required for the Wbb and Wcc 
samples by using events with two jets that pass the 
event selection described in Sec. Ullbut fail the 6-tagging 
requirements (zero-tag sample). The zero-tag sample has 
no overlap with the sample used to measure the single 
top quark cross section. During this study, we keep the 
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normalization of the lU+jets sample fixed to the value 
obtained by the iterative method described above and 
derive Ahf with the following equation: 



N (o) 



AO) 

'Wlp 



7 (0) 
*Whp> 



(2) 



whci'C N — jVdata — ^multijets — ^"non-Wjots, Ny/lp = 

iVwjj+iVwcj, andiVwhp = Mvcc+Mvbb- The superscript 
(0) indicates that the equation is written for the zero-tag 
sample defined above. The measured value of Ahf is 
consistent with one. Uncertainties on the assumed cross 
sections for single top quark, it, and Wcj production 
and the cross section ratio of Wcc to Wbb are taken into 
account. As expected, Xrf is most affected by varia- 
tions on the Wcj cross section and the Wcc to Wbb 
cross section ratio. An estimated uncertainty of 12% is 
assigned to the normalization of the Wcc and Wbb MC 
samples based on this study. 

We also consider other sources of systematic uncer- 
tainty from modeling both the background and signal. 
These uncertainties usually affect the normalization and, 
in some cases, also the shape of the distributions. The 
largest uncertainties arise from the jet energy scale (0.3- 
14.6)%, jet energy resolution (0.2-11.6)%, corrections to 
6-tagging efficiencies (6.6-21.2)%, and the correction for 
jet-flavor composition in lU+jets events 12%. There are 
also contributions due to limited statistics of the MC 
samples 6.0%, the measured luminosity 6.1%, and uncer- 
tainties on the trigger modeling 5.0%. 

Table Q] lists the numbers of expected and observed 
events for each process after event selection, including b- 
tagging. Figure [2] shows comparisons between data and 
simulation before and after applying 6-tagging. In the 
same figure, the normalization and differential spectra 
of the two dominant backgrounds are checked using the 
control samples dominated by VF+jets (e), and by it 
(f) events. These plots are indicative of the adequate 
background modeling attained for various sample condi- 
tions in the analysis. 

TABLE I: Numbers of expected and observed events in a data 
sample corresponding to 5.4 fb _1 of integrated luminosity, 
with uncertainties including both statistical and systematic 
components. The tb and tqb contributions are normalized to 
their SM expectations for a top quark mass of 172.5 GeV. 
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2 jets 
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4 jets 
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16 


44 
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13 
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3.5 


tqb 


140 
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13 


72 
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9.4 


26 
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6.4 


it 


433 


± 


87 


830 


± 


133 


860 
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163 


VK+jets 


3,560 


± 


354 


1,099 


± 


169 


284 
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76 


Z+jets & dibosons 


400 


± 


55 


142 


± 


41 


35 


± 


18 


Multijets 


277 


± 


34 


130 


± 


17 


43 


± 


5.2 



Sum of above sources 4,914 ± 558 2,317 ± 377 1,261 ± 272 
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FIG. 2: [color online] Comparisons between the data and 
the background model for (a) $t, (b) W boson transverse 
mass before 6-tagging, and (c) light quark jet pseudorapidity 
multiplied by lepton charge, after 6-tagging. Reconstructed 
top quark mass (d) after 6-tagging, (e) in a control sample 
dominated by VF+jets, and (f) in a control sample dominated 
by ti. The hatched bands show the ±1<t uncertainty on the 
background prediction for distributions obtained after 6-jet 
identification (c-f). The VF+jets contribution includes events 
from Z+jets and diboson sources. 



IV. MULTIVARIATE ANALYSES 

Since the expected single top quark contribution is 
smaller than the uncertainty on the background, wc 
use multivariate analysis (MVA) methods to extract 
the signal. The application of these methods to the 
measurement of the single top quark production cross 
section is described in Ref. [l7|. Three different MVA 
techniques are used in this analysis: (i) Bayesian 
neural networks (BNN) (fi) boosted decision trees 

(BDT) [4l[, and (iii) neuroevolution of augmented 
topologies (NEAT) [42[ . Each MVA method constructs a 
function that approximates the probability Pr(5'|x) that 
an event, characterized by the variables x, originates 
from the signal process, S = {tb, tqb, tb + tqb}. Therefore 
each method defines a discriminant D that can be used to 
constrain the uncertainties of the background in the low- 
discriminant region D w and extract a signal from an 
excess in the high-discriminant region Dwl, All three 
methods use the same data and model for background, 
performing the analyses separately on the six mutually 
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exclusive subsamples defined before. All three methods 
also consider the same sources of systematic uncertainty, 
and are trained using variables for discriminating signal 
from background chosen from a common set of well- 
modeled variables. These variables can be classified in 
five categories: single object kinematics, global event 
kinematics, jet reconstruction, top quark reconstruction, 
and angular correlations. The BNN uses four-vectors of 
the lepton and jets and a two-vector for $t to build 
the discriminant. The BNN performance is improved 
by adding variables containing the lepton charge and b- 
tagging information, resulting in 14, 18, and 22 variables 
for events with 2, 3, and 4 jets. The BDT ranks 
and selects the best fifty variables for all the analysis 
channels, while NEAT uses the TMVA HH implemen- 
tation of the "RulcFit" 44[ algorithm to select the best 



thirty variables in each channel. 

Each MVA method is trained separately for the two 
single top quark production channels: (i) for the tb 
discriminants, with tb considered signal and tqb treated 
as a part of the background, and (ii) for tqb discrimi- 
nants, with tqb considered signal and tb treated as a part 
of the background. 

Using ensembles of datasets containing contributions 
from background and SM signal, we infer that the corre- 
lation among the outputs of the individual MVA methods 
is ?s 70%. An increase in sensitivity can therefore be 
obtained by combining these methods to form a new 
discriminant @. To achieve the maximum sensitivity, a 
second BNN is used to construct a combined discriminant 
for each channel, for tb, tqb, and tb+tqb events, defined as 
B t b, B tqb and B tb+tqb . The B tb and B tqb discriminants 
take as inputs the three discriminant outputs of BDT, 
BNN, and NEAT, and they are trained by assuming tb 
or tqb as signals, respectively. The combined tb + tqb 
discriminant (B tb +tqb) takes as input the six discrim- 
inant outputs of BDT, BNN, and NEAT that are trained 
separately for the tb and the tqb signal. The training 
for Btb+tqb treats the combined tb+tqb contribution as 
signal with relative production rates predicted by SM. 
Figure [3] shows the outputs of the B tb , B tqb , and Btb+tqb 
discriminants, where good agreement is observed over 
the entire range. In these plots, the bins are sorted 
and merged ("ranked") as a function of the expected 
signal-to-background ratio (S:B) such that S:B increases 
monotonically within the range of the discriminant. For 
the tqb and tq+tqb discriminants, presence of signal is 
significant in the plots. For the tb discriminant, the signal 
presence is not as significant. 



V. MEASURING SIGNAL CROSS SECTIONS 

A. Bayesian approach 

We use a Bayesian approach [(| [TH, [l?} to extract 
the production cross sections. The method consists of 
forming a binned likelihood as a product of all six analysis 
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FIG. 3: [color online] Distributions of the (a) Btb, (c) Bt q b, 
and (e) B t b+tqb discriminants for the entire range [0-1] of the 
output. Distributions of the (b) Btb, (d) B tq b, and (f) Btb+tqb 
discriminants for the signal region [0.8-1]. The bins have 
been "ranked" by their expected signal-to-background ratio. 
The tb, tqb, and tb+tqb contributions are normalized to the 
measured cross sections in Table [TT] The hatched bands show 
the ±lcr uncertainty on the background prediction. 



channels (2, 3, or 4 jets with 1 or 2 6-tags) and bins 
using the full discriminant outputs. We assume a Poisson 
distribution for the number of events in each bin and 
uniform prior probabilities for non-negative values of the 
signal cross sections (tb, tqb and tb + tqb correspond- 
ingly). Systematic uncertainties and their correlations 
are taken into account by integrating over signal accep- 
tances, background yields, and integrated luminosity, 
assuming a Gaussian prior for each source of systematic 
uncertainty. A posterior probability density as a function 
of the single top quark cross section is constructed, with 
the position of the maximum defining the value of the 
cross section and the width of the distribution in the 
region that encompasses 68% of the entire area corre- 
sponding to the uncertainty (statistical and systematic 
components combined). The expected cross sections are 
obtained by setting the number of data events in each 
channel equal to the value given by the prediction of 
signal plus background. 
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B. Ensemble testing 

The methods used for extracting the cross sections 
are validated by studies performed using ensembles 
of pseudo-experiments that are generated taking into 
account all systematic uncertainties and their correla- 
tions. These ensembles of events are processed through 
each MVA method for each single top quark production 
mode and through the same analysis chain as used for 
the data. Five arbitrary signal cross sections (including 
the SM prediction) are used to calibrate the tb, tqb, and 
tb + tqb cross section extraction procedure. Means and 
standard deviations are determined by fitting Gaussian 
function to the distributions of extracted values of the 
measured cross sections in each ensemble. Figure 3] 
shows the resulting distributions and Gaussian fits for 
SM ensembles for tb, tqb, and tb+tqb processes. Straight- 
line fits of the extracted mean cross sections to the input 
values are shown in Fig.[5j where the shaded bands reflect 
the standard deviations of the extracted cross sections in 
each ensemble. 

The results of these pseudo-experiments show that the 
biases on the cross sections are negligible compared to the 
standard deviations of the extracted values. We therefore 
do not apply corrections to the measured values of the 
cross sections in data. 



section to vary substantially (as observed in the tb 
channel) or in a way that is not monotonic with the 
assumed top quark mass (as observed in the tqb channel). 



TABLE II: Expected and observed cross sections in pb for tb, 
tqb, and tb+tqb production. All results assume a top quark 
mass of 172.5 GeV. 
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C. tb, tqb and tb + tqb channel cross sections 

To measure the individual tb (tqb) production cross 
section, we construct a one-dimensional (ID) posterior 
probability density function with the tqb (tb) contri- 
bution normalized with a Gaussian prior centered on 
the predicted SM cross section and treated as a part of 
the background. This is implemented for each individual 
MVA method and also for their combination. To measure 
the total single top quark production cross section of 
tb+tqb, we construct a ID posterior probability density 
function assuming the production ratio of tb and tqb 
predicted by the SM. 

Figure [5] shows the expected and observed posterior 
density distributions for tb, tqb, and tb + tqb using the 
combined discriminants B t b, B tq b, and B t b+tqb, respec- 
tively. Table [TT| lists the expected and measured cross 
sections for the individual MVA analyses. All of the 
results are consistent with SM predictions, and the 
measured tb+tqb production cross section is the most 
precise current measurement, with a precision compa- 
rable to the world average ^M,. All results assume a top 
quark mass of 172.5 GeV and have a small correction 
for events with more than four jets based on the SM. 
The dependence of the measured cross section on m t 
is summarized in Table IIIII The assumed top quark 
mass affects the yield and differential properties for the 
signal acceptance and the modeling of ti events, which 
constitute the second largest background. The interplay 
between these two effects can cause the measured cross 



TABLE III: Dependence on mt of the measured cross sections 
in pb for tb, tqb, and tb+tqb production, using the combined 
discriminants for the assumed top quark masses. The 
predicted cross sections i in pb are also included in the table 
and labeled "SM". 
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VI. SIGNAL DOMINATED DISTRIBUTIONS 

Figurc[7]shows a comparison of the distributions of four 
kinematic variables with large discriminating power, for 
single top quark production in a data sample selected 
with S:B > 0.24 based on the B t b+tqb discriminant. 
Variables shown are: leading 6-tagged jets pr, W boson 
transverse mass, centrality, defined as the ratio of the 
scalar sum of the px of the jets to the scalar sum of the 
energy of the jets in the event, and reconstructed m*. 
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FIG. 4: Distribution and Gaussian fit of the measured cross section in a ensembie of pseudo-experiments with the same 
integrated luminosity as in data generated assuming the SM for (a) tb, (b) tqb, and (c) tb + tqb processes. 




FIG. 5: [color online] Mean (points) and standard deviation (shaded bands) of cross section as a function of the input cross 
section for the (a) tb, (b) tqb, and (c) tb+tqb single top quark processes from the ensemble studies of pseudo-experiments with 
the same integrated luminosity as in data. The continuous lines show the fits to the mean values where their uncertainties are 
smaller than the size of the points. The dotted lines represent the responses in the case of slope equal one and zero intercept. 
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FIG. 6: [color online] The expected (back) and observed (front) posterior probability densities for (a) tb, (b) tqb, and (c) tb-\ 
production. The shaded bands indicate the 68% C.L.s from the peak values. 
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The presence of the single top quark signal is needed to 
ensure a good description of the data. 

VII. \V tb \ MEASUREMENT 

The single top quark production cross section is 
directly proportional to the square of the CKM matrix 
element |Vtt>| 2 , enabling us to measure \V t b\ directly 
without any assumption on the number of quark families 
or the unitarity of the CKM matrix [l?} • We assume that 
SM sources for single top quark production and that top 
quarks decay exclusively to Wb. We also assume that the 
Wtb interaction is CP-conserving and of the V — A type, 
but maintain the possibility for an anomalous strength of 
the left-handed Wtb coupling (fi), which could rescale 
the single top quark cross section [45[. Therefore, we 
are measuring the strength of the V — A coupling, i.e., 
\Vtbfi\, which can be > 1. 



We form a Bayesian posterior \V t bf\\ 2 with a flat prior 
based on the B t b+tqb discriminant. Additional theoretical 
uncertainties are considered for the tb and tqb cross 
sections Q. Using the measured tb+tqb cross section, 
we obtain \V t bfi \ = 1.02+ 011 . If we restrict the prior 
to the SM region [0,1] and assume /f = 1, we extract 
a limit of \V tb \ > 0.79 at the 95% C.L. Figure [5] shows 
the posterior density functions for \Vtbf\\ 2 and for |Vth| 2 , 
assuming f[ = 1 and < |V t h| 2 < 1. 



VIII. SUMMARY 

In summary, we have measured the single top quark 
production cross section using 5.4 fb -1 of data collected 
by the DO experiment at the Fermilab Tevatron Collider. 
For m t ~ 172.5 GcV, we measure the cross sections for 
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Leading b-tagged jet p T [GeV] w Boson Transverse Mass [GeV] Centrality Top Quark Mass [GeV] 



FIG. 7: [color online] Distributions for data in the regions of large value for signal discrimination: (a) leading b-tagged jet pr, 
(b) W boson transverse mass, (c) centrality, denned as the ratio of the scalar sum of the pr of the jets to the scalar sum of the 
energy of the jets in the event, and (d) reconstructed nit- The contributions from signal have been normalized to the measured 
tb+tqb cross section. 
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FIG. 8: The posterior density functions for (a) \Vtbfi\ 2 and 
(b) |Vt(,| 2 . The shaded (dark shaded) band indicates regions 
of 68% (95%) C.L. relative to the peak values. 



tb and tqb production to be 

a{pp -+tb + X) = 0.68±g;|| pb 



a(pp -> tqb + X) = 2.86±g;|| pb 

assuming, respectively, tqb and tb production rates as 
predicted by the SM. The tqb cross section is consistent 
with the value cr(pp -> tqb + X) = 2.90 ± 0.59 pb 



measured in Ref. [27| . where we use the same dataset 
and discriminant but extract the cross section without 
any assumption on the tb production rate. The total 
cross section tb + tqb is found to be 



a(pp ->tb+ tqb + X) = 3.43 



"°- 73 nb 
-0.74 P D 



assuming the SM ratio between tb and tqb production. 
All measurements are consistent with the SM predictions 
for a top quark mass of 172.5 GeV. Finally, we derive a 
direct limit on the CKM matrix element |Vt&| > 0.79 at 
the 95% C.L. assuming a flat prior within < \Vtb\ 2 < 1- 
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